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Abstract 

We consider the anharmonic oscillator with an arbitrary-degree anharmonicity, a 
damping term and a forcing term, all coefficients being time-dependent: 

u" + gi(x)u' + g2{x)u + gz(x)u n + gi(x) = 0, n real. 

Its physical applications range from the atomic Thomas-Fermi model to Emden gas 
dynamics equilibria, the Duffing oscillator and numerous dynamical systems. The 
present work is an overview which includes and generalizes all previously known results 
of partial integrability of this oscillator. We give the most general two conditions on 
the coefficients under which a first integral of a particular type exists. A natural 
interpretation is given for the two conditions. We compare these two conditions with 
those provided by the Painleve analysis. 

1 Introduction 

The harmonic oscillator is the simplest approximation to a physical oscillator and, when 
perturbation terms are taken into account, the resulting anharmonic oscillator is governed 
by the nonlinear differential equation 

E = u" + gi u' + g 2 u + g 3 u n + g 4 = 0, n(n - l)g 3 ± 0, (1.1) 

where ' denotes the derivative with respect to the independent time or space variable x, 
gi(x) a damping factor, g 2 (x) a time-dependent frequency coefficient, gs(x) the simplest 
possible anharmonic term, g±(x) a forcing term. 

As to the anharmonicity exponent n, it can be either real if u(x) is real positive, which 
is the case for Lane-Emden [24J gas dynamics equilibria, rational, like n = 3/2 in the 
Thomas and Fermi |31[ [FT] atomic model, or more usually integer: —3 for the Ermakov 
[12j or Pinney [28] equation, 3 for the Duffing oscillator [10J. 

For generic values of the coefficients, this equation is equivalent to a third order au- 
tonomous dynamical system, which generically admits no closed form general solution. 
The purpose of this article is to review all the nongeneric situations for which there exist 
exact analytic results, such as a first integral or a closed form solution, either particular 
or general. This can only happen when the coefficients satisfy some constraints. 

The paper is organized as follows. In Section[2j we give a Lagrangian and a Hamiltonian 
formulation for any value of the coefficients (n,gi). This generalizes all the previous 
particular results, obtained for values of (n, g±, g 2 , g 3 , 54) equal to: 
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(5; 0, const,const, 0) [3], 
(5;2/x,0, 1,0) [261 Eq. (3.7)], 
(n; 51, 0,53,0) W\ E3 M , 
(n; 0, const, ax a , 0) [l], 
(n; pi, 0,1,0) [25], 

(n; 51, 52, 53,0) [33], [231 Section 6.74, vol. 1]. 

In Section [3l we provide two conditions on (n,gi) which are sufficient to ensure the 
existence of a first integral. 

In Section HI we give a natural interpretation of these two conditions. 

Finally, in section \5\ we perform the Painleve analysis of (jl.ip . Most of this work has 
already been done by Painleve and Gambier [19]. Indeed, the ordinary differential equation 
(ODE) (jl.ip belongs, at least for specific values of n and maybe after a change u i— > u N 
of the dependent variable u in case n is not an integer, to the class of second order ODEs 
which they studied and classified. However, as opposed to these classical authors, we do 
not request the full Painleve integrability of the ODE, only some partial integrability, 
and this requires some additional work. In particular, we compute the condition for the 
absence of any infinite movable branching, i.e. a multivaluedness which occurs at a location 
depending on the initial conditions. Such a condition, like for linear ODEs, arises from 
any integer value of the difference of the two Fuchs indices, whether positive or negative, 
and we check that this condition is a differential consequence of the two conditions for the 
existence of a particular first integral. This detailed Painleve analysis of equation (jl.ip 
happens to be an excellent example for several features of Painleve analysis which are 
most of the time overlooked. 

For convenience, we use the notation 



and the convention that function G\ implicitly contains an arbitrary multiplicative con- 
stant; letter K, with or without subscript, denotes an arbitrary constant. Function G\ 
frequently occurs, for the way to suppress term g\v! in (jl.ip is to perform the change of 
function u — > G 1 1 ^ 2 u. 

2 Lagrangian and Hamiltonian formulations 

For every value of (n, gi), including the logarithmic case n = — 1, the anharmonic oscillator 
can be put in Lagrangian form 





74 = Log 54 



(1.2) 




(2.1) 



or in Hamiltonian form 



q = p = 

op 



dH 
dq 



(2.2) 
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as shown by the explicit expressions for L, H, q,p 



u 



12 



2g 3 / u n du - g 2 u - 2g A u 
Jo 



u' 2 + 2# 3 / u n du + g 2 u 2 + 2g A u 



--hu, 



(2.3) 

(2.4) 
(2.5) 



L(u, u' , x) = G\ 

H(q,p,x) = Gi 

q = u, p = 2G\U + hu, 
in which h is an arbitrary gauge function of x. 

3 Particular first integral 

According to Noether theorem, one can find first integrals by looking at the infinitesimal 
symmetries of the Lagrangian. For a detailed review of this Lie symmetries approach to 
the anharmonic oscillator, the interested reader can refer to [14] . Since the dependence of 
ODE (jl.ip in u is rather simple, let us determine under which conditions on parameters 
(n,gi) there exists a particular first integral containing the same kind of terms than the 
Hamiltonian 

I = hu' 2 + h I" u n du + huu' + f A u 2 + hu + / 6 , (3.1) 
Jo 

in which the six functions fi of x are to be determined. 
Eliminating u" between /' and E, we obtain 

/' - (2/au' + hu)E = f 2 [ U u n du - g 3 f 3 u n+l + (fc - 2g 3 f 1 )u n u' + & (3.2) 

Jo 

+(/i + h - 251/1K 2 + + 2/4 - 51/3 - 2<? 2 /i V 
- 92h)u 2 + (h - 2g A h)u' + (f' b - g 4 f 3 )u. 

Out of the nine monomials u n du, n n+1 , u n u', u' 2 , u 2 , uu', u', u, 1, only eight are lin- 
early independent since n(n — 1) ^ 0, thus generating eight linear homogeneous differential 
equations in six unknowns, hence generically two conditions on (n,gi). Note that, even in 
the logarithmic case n = — 1, the first generated equation is f 2 ~ ( n + 1)53/3 = 0- Functions 
f 2 to / 6 are given by 



h 
h 



2<7 3 /i, 
25i/i -f{, 
■2 



3 1 

h = (9i +92- g'i)h - 2 5i/i + 2/1 > 

h = 254/1, 

/X 
gshdx, 

and function f\ must be a nonzero solution common to the three linear equations 
[_2(„ + l)g l93 + 2g' 3 ] h + {n + 3)g 3 f[ = 0, 

(-25152 + 25i5i - 5i + &)fi + (5i + ^2 - ^i)/i - \gxf'l + \f" = 0, 
(-25x54 + 25O/1 + 354/1 = 0. 



(3.3) 
(3.4) 

(3.5) 
(3.6) 

(3.7) 

(3.8) 
(3.9) 
(3.10) 
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Each above equation can be integrated once, 



pn+3 n-2n-2 „2 
53) 



K 2 



f?G^g 2 + 



{(gl-g[)h 



\gifi + \fD'hG? 



dx, 



(3.11) 
(3.12) 
(3.13) 



Whatever be (n, ft), the function can always be computed from ()3.9[) : depending on 
{ n -,9ij-, it is also given by 



n / -3 
54^0 
n = —3 



A 
/i 
h 



G 



2(n+l)/(n+3) -2/(n+3) 



.93 



.-,2/3 -2/3 



-5 3 v> 



general solution of 



1^3 / 



2 93 







(3.14) 
(3.15) 

(3.16) 



where the constants K\ and K$ have been absorbed in the definition of G\. The only case 
in which equation (|3.16|) needs to be considered is n = —3,54 = 0, and its solution can be 
found in [El Q2] [221 ^[14.33] [28] [U Eq. E12] [5]. 

Once fi is determined, f 2 and f§ are given by (|3.3|) . (|3.7|) . and /3,/4 by the three 
following expressions, corresponding to cases (|3.14p . (|3.15p . (|3.16p respectively, 

(25i+7 3 ) 



h 



h 



ra + 3 

^(251+74) 

-73 " T" 

, ~2(n-l)gf 
92 H 



(3.17) 



(n + 3)(2 5 ;+ 7 ^ 



n 



5)5173 + 27! 



(n + 3) 5 



52 + -5i 
1 



1 



5174 



274 



7^ 



(3.18) 



52^ 4 / 3 ^ 2 73 4 /3 /i 2/i 



Parameters (n,5j) must satisfy the conditions, polynomial in n, 51, 52, 73, 74, resulting 
from the elimination of /1 between the three linear equations (|3.8p - (|3.10p . There are two 
such conditions when 5354 is nonzero, and only one when it is zero. The simplest choice 
of these two conditions is (the labelling refers to the contributing 5i's): 



54 / : C134 = 2n5i - 37 3 + (n + 3)7 4 = 0, 

54 + : C 124 = 4 5 3 - I85152 - I851 + 275 2 + (65? - 36 52 + 275^)74 
-651 H ~ 4 7 4 3 + 95174 + 187474 " 974" = 



(3.19) 



(3.20) 



uniquely defined as, respectively, the condition independent of 52 and the one indepen- 
dent of (72,53). By elimination, one obtains the condition independent of 54 and the one 
independent of 51, 

C 123 = - 4 (251 + 73) 3 + (n + 3)[(n - 3) (-Agf - 2g'{ - 4g 2l ' 3 - 73") 
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34 + : C, 



234 



+ 2(n + 3)(n - 1)3132 + n(-83i3i - 63173 2 ) + (n + 3) 2 3 2 
+ (n - 9)(-2 5l 2 7 ^ - 3W3) " 3(n - 1)3173 + 67373] = 0, 

n% + n 2 (-^7 3 " 73-" + jpfri + ^ " 27 4 7 4 ' + iT) 



(3.21) 



3 /2 / 1 1 /3 2 / 1 \ / 
o7 3 7 4 + 7T^3 - n ( n - 1 )527 4 



^(- 2 - 3)7 3 74 2 + ^(n 2 -l)7f 



0. 



(3.22) 



For (n + 3)34 7^ 0, any two of the above four conditions are functionally independent. 
For n = —3, one has 27C123 ~~ 4Cf 34 = and independent conditions are C134 and C234. 
All above conditions admit an integrating factor, a natural consequence of the integrated 
forms (|3.1ip - (|3.13p . This is evident for C134; for each of the three others, it is sufficient 
to integrate it as a first order linear inhomogeneous ODE in 32 , 



34 ^ : if 134 
#4 / : if 124 
if 123 



2n — 3 „n+3 



sun -an 
^1 53 54 



32 - 



2 2 2 



gffl + ^3174 + g74 - g74 



G 



8/3-4/3 
i/4 ) 



[(n + 3) 2 92 - (n + 3)(25i + 73) - 2(n + 1)3? 



(n- 1)3173 +7?]" +3 G 2(w - 1) 



S3 



34 7^ : if : 



234 



in 



32 



+ 2)7 3 74 " 7 3 2 - (" + 1)7 4 2 + MH - 11) 



2n 2 



-4/3 4/n 
5 3 <?4 • 



(3.23) 
(3.24) 

(3.25) 
(3.26) 



In the Duffing case n = 3, condition C123 has already been given |15| . together with its 
integrated form if 123 [13] . 



4 Interpretation of the two conditions 

A very simple interpretation can be given for the two conditions. Indeed, the form of 
equation (|1.1|) is invariant under the simultaneous change of dependent and independent 
variables 



u(x) -> U(X) : u = a(x)U, X = £(x), 
where a and £ are two arbitrary gauge functions. The transformed ODE reads 



(4.1) 



u" + l - 

„n-l 



a 



31 



+- 



a' 



1 



a a 

32 +3i 1 

a a 



U 



t/2 



9zU n + —g 4 = 0. 



(4.2) 



Let us adjust the two functions a,£ so as to make two of the four new coefficients as 
simple as possible. One of the three possible ways is to cancel the damping term by the 
choice £' = aT 2 ^ 1 , which reduces ODE (H~2|) to 



U" + a G\ 



a' a" 

32 +3i 1 

a a 



U + a n+3 G 2 ig3 U n + a 6 G{g A = 0. 



3^2, 



(4.3) 
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Canceling the new g% coefficient amounts to solving the general linear second order 
ODE for a, which is possible (from the point of view of Painleve, adopted here) but does 
not lead to an explicit value of a. This reduced form is then 

U" + h 3 U n + h 4 = 0, (4.4) 

and this means that one can freely set g\ = g 2 = in (II. lj) without altering its global 
properties (existence of first integrals, Painleve property, etc). Instead of that, one can 
make constant either the reduced 53 coefficient iff (n + 3)<?3 7^ by choosing a n+3 = 
Gj 2 g^, or the reduced 54 coefficient iff #4 7^ by the choice a = G 1 2 ^ 3 g 4 1//3 (let us recall 
that Gi implicitly contains an arbitrary multiplicative constant). 
We are thus led to the reduced forms 



g A / : 0i 0, 54 1, 
92 1 



(4.5) 



92 



T01 



gSl + g9l74 + g74 



-3,54 = 



53 
91 

92 



54 
91 

92 



8/3 4/3 

W 9a > 

. „ -(n+3)/3 
> 53^1 54 ! 



0, 93 
[52 - 



1, 



(4.6) 



; (H +73) + 



: (-2(n + !)</? 



n + 3 v J1 ' ' (ra + 3) 2 

/- -n '1 «M/-f2(n-l)/(n+3) -4/(n+3) 
- 1)5173 + 73 )] G i 5 3 7 

^271/(71+3) -3/(n+3) 
54^1 93 ) 

0, 53 >->• 93GI, 

a' a" 

92 +91 1 

a a 



a A G\ 1 ^ 0. 



(4.7) 



Then the interpretation is obvious: any reduced coefficient distinct from or 1 is 
the r.h.s. of one of the integrated conditions (|3.23j) - ()3.26p . Conversely, any integrated 
condition is one of the remaining coefficients when two coefficients have been made constant 
by a choice of gauge. For instance, -K234 is the reduced 52 coefficient associated to reduced 
coefficients 53 and 54 equal to unity. 

This can also be seen in a more elementary way. In a gauge (a,£) such that g\ = 
0,g' s g' 4 = 0, an expression for the first integral is 



91 = 0, g 3 g 4 = : I 
and, from the relation 



u' 2 + 2 ff3 



u n du + g2U + 2<74ii, 



I'q - 2u'E = 2g' 3 / u n du + g' 2 u 2 + 2g' A u, 
Jo 



(4.8) 



(4.9) 



one deduces that the two other coefficients g 2 and g 3 or 54 must be constant. 

The Hamiltonian (|2.4p is a first integral if and only if g\ = and all other g^s are 
constant. 
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5 Painleve analysis 

Painleve set up the problem of finding nonlinear differential equations able to define func- 
tions, just like the first order elliptic equation 

u' 2 = 4-u 3 — giu — (?3, (§2,93) complex constants, (5-1) 

defines the elliptic function of Weierstrass p(x, g%, gs), a doubly periodic function which 
includes as particular cases the well known trigonometric and hyperbolic functions. For a 
tutorial introduction, see the books [21\ 17]. 

A by-product of this quest for new functions has been the construction of exhaustive 
lists of nonlinear differential equations, the general solution of which can be made single- 
valued (in more technical terms, without movable critical singularities, this is the so-called 
Painleve property (PP)), which implies that their general solution is known in closed form. 
In particular, the list of second order first degree algebraic equations, i.e. 

u" = F(u',u,x), (5.2) 

with F rational in u', algebraic in u, analytic in x, which possess the PP has been estab- 
lished by Painleve and Gambier [19J. 

These classical results apply to our problem only for those values of n for which 
Eq. (ll.ip . maybe after a monomial change of the dependent variable u = U k ,k G 1Z, 
belongs to the class (I5.2D . These values, which include at least all the integers, are deter- 
mined below. Then, the way those classical results can be applied is twofold. 

1. Require the PP for our equation or its transform under u = U . 

2. Restricting to the values of (n, gi) for which the first integral (|3.ip exists, check that 
the two conditions for the existence of this first integral imply the identical satisfac- 
tion of the necessary condition that Eq. (|5.2p have no movable logarithmic branch 
points. Indeed, this is a classical result of Poincare that the movable singularities 
(i.e. those which depend on the initial conditions) of first order algebraic ODEs 
can only be algebraic, i.e. u ~ uq(x — xq) p , and never logarithmic, i.e. with some 
Log(x — xq) term. Let us do that without too many technical considerations. 

The above mentioned necessary condition that Eq. (|5.2p have no movable logarith- 
mic branch points can only be computed after performing the following steps (for the 
unabribged procedure, see [U section 6.6]). 

Step 1. For each family of movable singularities 
+00 

j=0 

determine the leading behaviour (p,uq). This is achieved by balancing the highest deriva- 
tive u" with a nonlinear term. Therefore, there exist two leading behaviours, denoted 
"family 53" (balancing of u" and g3U n ) and "family 54" (balancing of u" and 54) 

1 

(93) ■ P = 7, u = 

n — 1 

(g 4 ) : p = 2, n = ~l;94, 54 + 0. (5.5) 



n + 1 
(n — \y 



n — 1 



(5.4) 
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Step 2. For each family, compute the Fuchs indices, i.e. the roots i of the indicial 
equation of the linear equation obtained by linearizing (jl.ip near its leading behaviour 
u ~ uqx p , an d require every Fuchs index to be integer. This linearized equation is 

(53) : (^+n< 73 ta 2 r- 1 )« = 0, (5.6) 

(54) : ( d ^ + °) w = ' ( 5 - 7 ) 
and the Fuchs indices are obtained by requiring the solution v = i>oX p+ *(l + 0(x)), 

(93) : (i + 1) (i ~ 2 - -^—) = 0, (5.8) 



n — 1 / 

( 9i ) : (i + l)(i + 2) = 0. (5.9) 

The diophantine condition that z = 2 + 4/(n — l)be integer has a countable number of 
solutions since we have not yet put restrictions on n. 

Step 3. For each family, compute all the necessary conditions for the absence of mov- 
able logarithms (in short, no-log conditions), which might occur when one computes the 
successive coefficients Uj of (|5.3p . One can check that the family g± can never generate 
such no- log conditions. These conditions need not be computed on the original equation 
(|l.ll) . they can be computed on any algebraic transform if this proves more convenient 
(indeed, movable logarithms are not affected by an algebraic transform on u), such as 

u = U k : kUU" + k(k - 1)U' 2 + k gi UU' + g 2 U 2 

+ g 3 U 2+(n - 1)k + g 4 U 2 - k = 0. (5.10) 

The transformed powers p are p 3 = — 2/((n — l)k),p4 = 2/k, and the Fuchs indices are 
unchanged. 

The computation of the no-log conditions is impossible unless there exists a k making 
all the powers of U in (|5.10p at least rational. In order to avoid the technical complications 
of dealing with rational values of the leading exponent p, we restrict to those values of n 
for which there exists a k making 2 + (n — l)k and, if 54 is nonzero, 2 — k integer. The 
useful transforms are 

u = v 2 '^ : -^—vv" - 2^—^v' 2 + -^— 9l vv> + g 2 v 2 
n — 1 (n — \) z n — Y 

+ 9 3V 4 + g4V 2 - 2/{n ~ 1) =0, (5.11) 

1 /(n — ]) 1 // ^ 2 ,9 1 1 o 

u = w /y ' : ww — -ttW H g\ww + g 2 w 

n — 1 (n — iy n — 1 

+W+W" 1/(n ~ 1) =0, (5.12) 
u = V~ 2 : -2VV" + 6V' 2 - 2 gi VV' + g 2 V 2 

+g 3 v 4 - 2n + g4V 4 = o, (5.13) 

u = W~ 1 : -WW" + 2W' 2 - giWW' + g 2 W 2 

+9sW 3 ~ n + 9i W 3 = 0, (5.14) 
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which are polynomial if and only if 



(EH]) 



54 = or (<?4 ^ and 



ez) 



dsn 



94, 

2n € Z, 
n£ Z. 



n — 1 

or (g 4 ^ and € -Z), 



n 



1 



(5.15) 

(5.16) 

(5.17) 
(5.18) 



The original ODE (fLTj) is identical to (j5TT|) for n = 3 and to ([5TT2]) for n = 2. 

To summarize, let us compute the no- log condition Qi = on the ODE for v (|5.1ip . 
Unfortunately, one does not know how to obtain the dependence of Qi on n, since re must 
first be given a numerical value before Qi is computed; this makes uneasy the comparison 
with conditions (|3.19p ^ (|3.20p . which depend on re. 

To fix the ideas, a list of useful values of (n,i) is displayed in Tabled! 



Table 1. Values of (i,n) for i integer £ [—4, 10]. 



4 

2 + 7 

re — 1 


-4 


-3 


-2 


-1 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


n 


1/3 


1/5 





-1/3 


-1 


-3 


oo 


5 


3 


7/3 


2 


9/5 


5/3 


11/7 


3/2 



The computation of Qi for positive values of i is classical |191l2]. Denoting for shortness 
C\ = Ci23,C2 = C134, one finds the following expressions Qi for the indicated values of 

(n,g4), 

(-3,0) : Qi = Ci, (5.19) 
(5,0) :Q 3 = Ci, (5.20) 

(3, 54 ) :Q 4 = i^I^r 1 / 2 !^- 5i)Ci-C(] + ic 2) (5.21) 

(1,0) : Q 5 = (-5? + -92 + .<?! + -5173 + -l'i)Ci + (-91 + -Vs)^ + -Cf, (5.22) 

(2, 54 ) : Q 6 = .d + .C? + .C( + .C'{ + .C(" + .C 2 + .C^, (5.23) 

(|,0):Q 7 = .Ci + ... + .cf\ (5.24) 

(|,0):Q 8 = .d + .-. + .C^. (5.25) 

where dots stand for rational numbers when i = 5 and polynomials of g\, g2, 73, 74 when 
i > 5. Similar relations have been checked for i = 9 and i = 10 (Thomas-Fermi case) 
but are not reproduced here. Condition Q4 = contains a ± sign arising from the two 
possible choices for i>o and is equivalent to the two conditions 

(73 - 5i)Ci23 - C[ 23 = 0, C 134 = 0. (5.26) 

We therefore check the property that each Qj is indeed a differential consequence of 
the two conditions C123 = 0, C134 = for the existence of a first integral (|3.ip 

Vz G AA,V 5i : (Ci = 0,C 2 = 0) => (Qi = 0). (5.27) 
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For negative [THl [6] values of the Fuchs index i, the results [6] are the following: the 
family never generates any no- log condition, and, for the family g^, a no- log condition 
arises from the Fuchs index —1, and this condition is a differential consequence of condi- 
tions (|3,19p - (|3.20p . at least for the examples handled (n, r, #4) = (1/5, —3,0), (1/3, —4,0). 
This is also an experimental verification of 



and this relation cannot be reversed, as proven by Painleve and Gambier. For instance, 
in the case of the Duffing oscillator (n,i,g&) = (3,4,54), condition Q4 = implies the 
reducibility of v to the second Painleve transcendent whereas the stronger conditions 
C\ = 0, C2 = imply the reducibility of v to an elliptic function. 

Remark. When one includes the contribution of the Schwarzian in the definition of the 
gradient of the expansion variable Xi as done in the invariant Painleve analysis [3], 



all the computed no-log conditions Qi = 0, equations (|5.19j) - (15.25p . are independent of 
this Schwarzian S, as opposed e.g. to the Lorenz model [9]. This certainly indicates some 
hierarchy between the level of nonintegrability of these two dynamical systems. 

Remark. For some small values of there is equivalence between the no-log condi- 
tion and (|3.19p - (|3.20p . This nongeneric situation occurs only for the following values of 

(M,S4), 

(—3,1,0), i.e. the Ermakov-Pinney equation [121 [25], 

(5,3,0), i.e. an equation considered by Lane and Emden [2H[TT], Chandrasekhar [3] and 
Logan [26, p. 52], 

(1/5, —3,0), an equation which could deserve more study. 

6 Conclusion 

This work generalizes all previous results on the partial integrability of the anharmonic 
oscillator. It gives a natural interpretation of the two conditions for the existence of a 
particular first integral, in terms of reduced coefficients. Finally, this system is an excellent 
example to study several features of Painleve analysis. 

A good, recent bibliography can be found in Ref. |20j . 
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Vi 6 Z, V gi : (Ci = 0, C 2 = 0) => (Q_ a = 0) 



(5.28) 



x' = i + fx 2 



(5.29) 
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